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Abstract 



We calculate the next-to-leading order short distance QCD corrections to the coefficient 
rji of the effective AS = 2 hamiltonian in the standard model. This part dominates the 
short distance contribution {Amx)^^ to the Kl - Ks mass difference. The next-to-leading 
order result enhances rji and (Am^^)^^ by 20% compared to the leading order estimate. 
Taking 0.200 GeV < A^jg < 0.350 GeV and 1.35 GeV < mc(mj < 1.45 GeV we obtain 
0.922 < r]f^^ < 1.419 compared to 0.834 < t]^^ < 1.138. For Bk = 0.7 this corresponds 
to 48 - 75 % of the experimentally observed mass difference. The inclusion of next-to- 
leading order corrections to rji reduces considerably the theoretical uncertainty related to 
the choice of renormalization scales. 
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1 Introduction 



Over the last three decades the study of the K-meson system has contributed enormously to 
our insight into the fundamental principles of nature, just to mention the discovery of CP- 
violation by Christensen et al. [|l| and the conjecture of the existence of the charm quark by 
Glashow, Iliopoulos and Maiani . The latter was derived from the suppression of flavour- 
changing neutral currents, which is responsible for the smallness of K'^ — if°-mixing. In 
the following years this AS* = 2 weak interaction process has been used to investigate 
the indirect CP-violation apparent in the Kl — Ks mass difference and to constrain the 
parameters of the Standard Model. The first highlight was clearly the prediction of the 
mass of the charm quark prior to its discovery by Gaillard and Lee ||^. 

The theoretical determination of the effective low energy AS = 2-hamiltonian, however, 
is not satisfactory yet, because the indispensable incorporation of strong interaction effects 
is difficult due to the poorly known infrared structure of quantum chromodynamics. The 
theorist's tool to tackle the problem is the use of Wilson's operator product expansion 
factorizing the Feynman amplitudes into a long distance part, which has to be calculated by 
non-perturbative methods, and a short distance Wilson coefficient, to which perturbative 
methods are applicable. The largely separated mass scales in the problem generate large 
logarithms, whose resummation is mandatory. This is achieved by renormalization group 
methods applied to a cascade of effective field theories, in which the heavy degrees of 
freedom are successively removed. 

Pioneering work in the short distance calculation of the AS = 2-hamiltonian was done 
by Vysotskij and by Oilman and Wise |0, who have used the leading logarithmic (LL) 
approximation. This method leaves certain questions, which will be summarized below, 
unanswered, and therefore many attempts have been made to include subleading effects. 

To be specific, consider the hamiltonian: 
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where is the Fermi constant, Xj = VjdVj's denotes the CKM-factors, and Qll is the 
local four quark operator (see fig. |^ on p. [TO ) 



Qll = (sj7/i(l - 75)c?i)(sfc7^(l - 75)4) = {sd)v-A{sd)v-A (2) 

with j and k being colour indices. 

In eq. (|l]) the OIM mechanism A„ + Ac + At = has been used to eliminate A^. S{x, y) 
and S{x) denote the Inami-Lim functions resulting from the evaluation of the box 
diagrams of fig. |l|. They must be understood to be expanded to lowest nonvanishing order 
in mil Mir in (1). 

The rjiS contain the QCD corrections with their dependence on the renormalization 
scale IX explicitly factored out in the coefficient &(/i). In the LL approximation 6(/i) reads 
in terms of the three-flavour strong interaction constant: h{ji) = [a3(yu)]~^^^. In the absence 
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Figure 1: The lowest order box diagrams contributing to H^^ ^. The zigzag hues stand 
for W-bosons or ficticious Higgs particles. 

of QCD corrections rjih{n) = 1. Since physical quantities must not depend on fi, b{fi) has 
to cancel the corresponding dependence in the matrix element of Qll, which reads in the 
case of K'^ — iiT^-mixing 

(k^IQlMIk'^) = IbmiWk- (3) 

Here niK and fx are the mass and decay constant of the i^'-meson, and the 'bag'-parameter 
Bxifi) combines with to the renormalization group invariant 

Bk = BK{fi)b{fi). 

In the following we mention the main reasons for going beyond the LL approximation 
in the calculation of rji. 

• The running charm quark mass rric enters H^^^^ at the scale /ic — rric, at which the 
effective four-quark theory is matched to an effective three-quark theory. The LL 
result depends considerably on fic- 

• Similarly rji depends on the scale nw ~ Myy, at which the W-boson is removed from 
the calculation. 

• The definition of the QCD scale parameter A requires a next-to-leading order cal- 
culation, and therefore one can at most quahtatively discuss the dependence on A of 
physical quantities calculated in the LL approximation. 

• After all the subleading terms can be sizeable. The LL hamiltonian only explains 
about 60 % of the observed — Ks mass difference. 

Buras, Jamin and Weisz ^ have improved the result of by including next-to-leading 
order (NLO) effects into 772 for if ^•5=2 ^^d H^^^^, which, being sensitive to the unknown 
top-quark mass, is relevant for indirect CP violation in K-decays and — i?*^-mixing. 
They have thoroughly verified that their Wilson coefficient depends neither on the gauge of 
the gluon propagator nor on the infrared structure of the factorized Feynman amplitude and 
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proved previous attempts to incorporate NLO effects into H^^ ^, H^^ ^ to be incorrect. 
Since the present worlc extends theirs to the coefficient rji, we closely stick to the notation 
in 0. 

An earlier attempted NLO calculation of rji has already been commented on in 
1^. The authors of ||^ have incorrectly extracted the short distance part of the transition 
amplitude and have thereby found a drastic decrease in rji rather than the sizeable increase 
found by us. 

The coefficient rji plays the key role for the Kl — Ks mass difference, while r]2 and 
?73 contribute only less than 10%. Since the LL calculation has yielded a too small mass 
difference, the correction in 0] even worsens the LL result. 

The nonperturbative parameter has been calculated using lattice gauge theories 
0, l/N expansion Q and QCD sum rules leading to Bk = 0.7 ± 0.2. The QCD 



hadron duality approach favours lower values B^ = 0.4±0.1 |13|. For a large calculation 
of both short and long distance contributions to the — A's-mass difference see M . 

Our paper is organized as follows: In section 2 we will set up our notations and present 
the general outline of the calculations. Section 3 briefly summarizes the LL calculation of 
changing the assumption nit <^ into the present day relation mt>M]v. Section 4 
is devoted to our NLO calculation followed by section 5 containing the numerical analysis 
of our results. We close the paper by summarizing our main findings. 



2 Preliminaries and Conventions 

Throughout this paper we work in the MS-scheme using an arbitrary _Rg-gauge for the 
gluon propagator and the 't Hooft-Feynman gauge for the W-propagator. Infrared (mass) 
singularities are regulated by quark masses. Encouraged by and |16| we use an anti- 
commuting 75 (NDR scheme). 

The algorithm to perform NLO QCD corrections to weak processes has been outlined 
ISj p!7| , 1^, so that we just have to summarize briefly the necessary steps with emphasis 
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on the new features in our calculation. 

The part of H^^^^ in (|ID relevant for rji reads 

H'^if^) = ^M^A^r/l(/iH.,/x.,/ic)■^(47^^)■%)■Q^^(/^)' (4) 

where we have stressed that the calculation of rji involves the scales fiw — Mw, l^b — rnb 
and /ic — rric, at which the W-boson, the bottom and the charm quark are removed from 
the theory. A physical quantity cannot depend on these scales and any residual dependence 
of H'^ on them is due to the truncation of the perturbation series. To judge the 'theoretical 
error' associated with this we have varied these scales instead of setting them exactly equal 
to the values of the corresponding current masses. We will see that the NLO calculation 
reduces the dependence of the LL result on fj,^, /i^ and fic- 
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2.1 The Hamiltonian at the Scale fiw — 



Consider first the Green's function witli tlie insertion of four weak currents calculated in 
the full standard model. The lowest order diagrams contributing to it are shown in fig. |I], 
the 0(a)-corrections are displayed in fig. 0. 





Figure 2: The classes of diagrams contributing to the 0(a;)-correction to {H'^)\ the 
remaining diagrams are obtained by left-right and up-down reflections. The curly lines 
denote gluons. Also QCD counterterm diagrams have to be included. Diagram Fg evaluates 
to for zero external momenta. 

After multiplication by i they constitute the matrix element of between quark states 
at the scale nw — M^: 

mfiw)) = {HM)^'^ + ^^{HMy'^ +Oia'). (5) 

In view of the fact that the desired short distance effects do not depend on the treatment of 
the external states, one can choose free quark states with zero external momenta through- 
out the whole calculation. The diagrams of fig. |1| with internal quark masses mj, rrij involve 
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S{mf/M^,myM'^) @, 0, which enter the coefficient of Tji in the GIM-combination (with 
0) 



2 2 

~, mt mt 



)-2S{ 



2 2 

~, mz m 



l) + S{ 



2 2 

~, mf, mf, 



S{ 



to give 0: 



(6) 



(7) 



The diagrams of fig. ^ have been calculated in giving after inclusion of the GIM mech- 
anism: 



(1) 



^ F \ 2 2 
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Only the lowest order in ml/M^ <C 1 has to be kept in (|^) and 
In (H) new operators have emerged: 



N 



2N ' 
'N'^ + N -1 
2iV 



1 



/jV^ + _ 1- 1 



(9) 



where L = 1 — 75,/? = I + 75, A^is the number of colours, and 1 and 1 denote colour singlet 
and antisinglet, i.e. (L ® i?) ■ 1 stands for Si{l — 'y^jdj ■ + 75)^4. We have written JJ 
in a manifestly Fierz self-conjugate way. The Dirac structures of Qll (see (|^)) and T are 
also Fierz self-conjugate. 
The functions in (P) are: 
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+ 



N 



2N 



61og^ + n2-2 

-ll + -7r2 + 121og-f + 31og 
3 Mr 



ml log(mf / /i^) - ml \og{my /i^ 



m^ml 



61og^ 



hu 



(-3-0 

3 -|- ^ m^mg , 



ml \og{ml/ /i^) - log(my /i^ 



mt, 



log- 



2 ' 



(10) 



5 



where Cp = (N'^ - l)/i2N). 

Let us now discuss the structure of {H^ny/)) in equations (|^) to (|l^): Notice first 
that there is no large logarithm in (|^, the logarithms of the individual box diagrams are 



cancelled by the GIM mechanism. In hiL in ([T0|) we have arranged the logarithms such 
that one can easily distinguish those which are small for /i ~ irtc from the large logarithm 
log (m^/M^). In the LL approximation the terms of the form [a;/(47r) ■ log(m^/M^)]"' , n = 
0, 1, 2, . . . are summed to all orders in perturbation theory. 

Next one observes that {H {nw))^^^ involves the additional operators U and T, which 
are artefacts of the use of small quark masses and m<j to regulate the mass singularities. 
Since the desired Wilson coefficient must neither depend on the gauge parameter ^ nor on 
the infrared regulators, hu, hx and the gauge part of Hll must disappear from it, thereby 
providing a check of the factorization described in section |2.2| . Yet in ( |T0| ) one can see that 



these terms do not depend on M\y, which sets the scale for the short distance physics. 

A dimensional regularization of the mass singularities would have resulted in hu = 
hx = 0. 

2.2 Factorization of the Hamiltonian at the Scale fiw 

In order to use the renormalization group (RG) evolution to sum the large logarithm 
log {m'^/M^r) to all orders in perturbation theory we seek a factorization of the form 

oc ^CjCj (11) 

j 

with composite operators Oj describing the weak interaction in an effective field theory 
without W-boson and top quark. The short distance physics is contained in the Wilson 
coefficients Cj, to which renormalization group methods will be applied. 

Consider first the time ordered product of two weak currents Jfj.{x) folded with the 
W-propagator T>'^'^, which mediates a AS* = 1 transition. Its operator product expansion 
is well known to order a |TB|, |T3|, |T^: 



-ijdSTj;^{x)Jt\yW{x-y) = C+g?(x) + C_gi'^(x) + . . . (12) 
Here j, k are fiavour indices standing for u or c. 

Q'i = {sj)v-A{kd)v-A-\(l±l) 

are dimension-six four-quark current- current operators, the meaning of 1 and 1 has been 
explained after eq. @). The dots in (|I^ ) denote penguin operators, which contribute to 773 



rather than 771, and operators with dimension higher than six. 

Now the Green's function described by the diagrams of fig. |l| and ^ contains two in- 
sertions of ([I2D and thus the effective theory involves matrix elements with two insertions 
of Qj^ or Q_. In order to find the complete set of operators needed in (|Tl|) we have to 
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look at the forthcoming renormahzation group analysis now: The reader may have noticed 
that we have already put the current-current operators into the basis {Q^.QJ), which is 
diagonal with respect to the renormahzation group evolution. Further the GIM mechanism 
prevents the mixing under renormahzation with penguin operators 
For brevity we define the bilocal operators 



m 



-QT{x)Qf{y)-QT{x)Qf{y) 



(13) 



The diagrams contributing to the matrix element {Oij) are depicted in fig. ^ and ^ the 
cross denotes the insertion of Qi or Qj. 




Figure 3: Diagram Dq in the effective four-quark theory constituting the matrix element 
of Oij to order The cross denotes the insertion of an effective AS = 1 operator. 

Now one expects the diagrams of fig. |^ and ^ to be divergent enforcing an extra renor- 
mahzation of Oij proportional to a local four-quark operator. The GIM mechanism, how- 
ever, cancels the divergences. Clearly in fig. ^ one-loop QCD counterterms must be in- 
cluded first. So we end up with 

H%f,) = ^Xl Y: C.jif^)0^j{f^). (14) 

^ i,j=+- 

These features have been discovered by Witten ||l9l a long time before we have carried 
out this explicit calculation. He has considered four-quark QCD and states: 

H%x) = ^Jd'yT 

where we have denoted the part of the AS = 1 hamiltonian containing the first two 
quark generations by H^'^^^^. For our case this implies that we can simply obtain Cij 
by multiplying the AS' = 1 Wilson-coefficients Cj and Cj calculated in |T^, where the 
AS = 1 hamiltonian is normalized as 

V ^ i=+,- kl 



H'^^^=\x)H'''^^=\y) , 
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Figure 4: The classes of diagrams in the effective four-quark theory contributing to 0{a). 
The other members of a given class are obtained by left-right and up-down reflections. 
Also QCD counterterms have to be included. Diagram Dg = for zero external momenta. 

This multiplicativity of Wilson coefficients has also been proven for a scalar theory by Lee 
| 20| . Since we anyway need the finite parts of the diagrams of fig. ^ and ^ later in the 
calculation (section we have compared J2Ci{fJ,w)Cj{^w){C^ij{fJ'w)) with {H'^i^ny/)) as 
given in eq. (H) and (H) and have indeed verified 



(15) 



to order a. 



2.3 Evolution of the Hamiltonian from fiw to /i^ ~ rric 

In the next step we want to evolve Cij{fi) from fj,\Y down to fic using the renormalization 
group formalism. This procedure is well known ^ [l^ and we briefly list its main 
ingredients: 

The running of the coupling constant g {a = g^/ (47r)) is determined by the beta function 



Pig) 



-A 



(/) 
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PI 
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0{g' 



with / being the number of flavours and 



(16) 



(17) 
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a{ii) reads in terms of the f-flavour QCD scale parameter A/: 



a(/i) 
47r 



/3(^hog(/.VA2) (/3(^))^log2(;,7A2) 



+ 



'logV] 



A3 and A5 are obtained from A4 = Aj^ by imposing continuity on a. 

The anomalous dimension 'jij of Oij are readily obtained from those of Q± via their 
definition in terms of the renormalization constants Z±. Since no matrix renormalization 
is required for Q± and the bilocal operators Oij do not mix with local operators, also Oij 
renormalizes multiplicatively, so that 



7ij 



The coefficients 7^°'' and 7^'^''^-''' of the anomalous dimension of Q±, 



(ZiZj) ^ fi— [ZiZ. 



li + lj, hj = +,-- 



1±{9) 



(0) 9 

7± — 
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read in the NDR scheme 15 



(0) 

7± 



(i)(/) 
±6 — — — , 7± 



We define for i,j = +, — : 



'(/) 
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(0) 



d. 



if) 



2A 



(/) 



7U 



if) 



if) , yif) 



zY' + z 



3 ' 



(19) 



(20) 



so that the solution to the renormalization group equation for Cij{fi) reads below the 
b-quark threshold fib in NLO: 



^(5) 



-,(4) 



a(/ifc) / \ a{iJ,) J 
^ ajfiw) - «(/^&) ^(5) ^ a(/ib) - a(/i) ^(4) 



An 



An 



Hence the RG improved Wilson coefficient Cij also equals Ci ■ Cj. 



(21) 
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2.4 The Hamiltonian below fi^ 



Below the scale /Xc — fTic the c-quark field is no more a dynamic degree of freedom and we 
have to express H^{^) in an effective three-quark theory. Here we are left with the single 
local four-quark operator Qll defined in eq. (0). 
To obtain the Wilson coefficient C in 

^^(/^) = ^^lC{f,)QLM (22) 
one must first solve the matching condition between the four- and three-quark theory: 

^Xl Y: QMiO.M) = ^,XlC{f^^){Q,M). (23) 

At this stage the matrix elements {Oij) must be calculated. (Qll) has already been 
obtained to order a in |^. The corresponding diagrams are depicted in fig. ^ and |^. 




Figure 5: The diagram for the matrix element of Qll in the effective three-quark theory 
to order a^. The cross denotes the insertion of the effective AS* = 2 operator. 




Figure 6: Classes of diagrams in the effective three-quark theory contributing to {Qll) 
to order a. The other members of a given class are obtained by left-right and up-down 
reflections. QCD counterterms have to be included. 

Since the LL analysis sums up the terms [a ■ large log]" for = 0, 1, 2, . . . , to all orders 
in perturbation theory, the LL matching has to be made to order a^. Consequently in 
NLO the terms a [a ■ large log]"', n = 0, 1, 2, . . . , are summed, and one has to match the 
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effective theories to order a^. For this the finite parts of the diagrams of fig. ^ must be 
calculated. 

This is in contrast to the recent NLO analysis of —>■ ir'^i'V, jJ^^ ^-^^ by Buchalla 

and Buras which also involves the renormalization of bifocal operators. There the one- 
loop Green's function already contains a large logarithm and therefore the NLO matching 
can be made at the one-loop level. 

The evolution of C (yu) in ( p2|) to scales below is dictated by the anomalous dimension 
of Qll, which equals 7+ in (|19|) |[T5|] : 

^{3) 

C(.) = C(,.) (^) * (1 + -AtAzJ^^zf + 0(«^)) . (24) 
The residual /x-dependence in (|^) is absorbed into defined in (P: 

%) = (l-^zf + 0(«^)). (25) 

We close this section by noting that C(/i) has via (|23|) absorbed all dependence on the 
charm quark mass present in {Oij{fic)), because the matrix element {Qhhi.lJ'c)) is evaluated 
in an effective three-quark theory containing no information on irtc- Consequently C(/i) 
contains the running charm quark mass in the form mdfJ'c)- 

To discuss the dependence of our result on the scale /ic in section we need the 
anomalous mass dimension 7^ , which is given by 



2 / ^2 



2 



with 



Hence with 



7i°) = GC„ ^(^^(f^ = ( + - H/ ) . (26) 



3 3 



dif) = ^ and Z^^ 



the running mass reads 



mc{lic) = mdmc) — r 1 + Z^i> . 27 

\a{mc) J \ An J 
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3 The Leading Log Analysis 



In this section we review the LL analysis of using the formahsm of section ^ and taking 
into account that the top quark is heavy. This analysis has already been performed by 
Flynn [Q. A good approximation neglecting quark thresholds has been given by Datta, 
Frohlich and Paschos ||23|. We discuss the scale dependences inherent to the LL calculation. 



The relevant RG formulae are obtained from those in section ^ by setting [3i 
7m '"■'^^ to zero. 



(/) 



7£^'''^'' and 



3.1 The Hamiltonian between the Scales and 

At first we want to factorize the H'^{iJ,) for jJ^c fJ' ^ fJ'W in order to separate short and 



long distance physics as described in sections ^]2| and [27^ : 



ry2 



8|) simply equals {W^i^nw)) 



(0) 



f28) 



given in (J^) . 



At the matching scale the left hand side of 

The AS" = 1 Wilson coefficients C±{^\y) are equal to 1 in the LL approximation [jl5 
so that we get Cij^Hw) = 1 by means of ([T5|). 
In section we will need (see fig. ^ 



'I] 



(29) 



with 



N + 3 



-A^ + 1 



Using ( [29| ) and (|^ one easily checks (p^). 

{H^{ixy/)) contains the Inami-Lim function S*, which has been expanded to lowest order 
in m'l/M'yY in (^). To include higher order terms, (m^/M^)^ ... in the factorization ([25|), 
one must take into account operators with dimension higher than six in (|12]). 

Next we seek the resummation of the large logarithm a log (m^/M^) present in the 
higher order terms of H'^. The factorization ([28| ) splits log (m^/M^) into log (/i^/M^) + 
log(m^//i^). While the latter resides in the matrix elements, the former is contained in 
the Wilson coefficients and is summed to all orders by the renormalization group evolution 
from /ivi/ to /ic- 

The LL running of the Wilson coefficients down to is given by 



a(/ic 



^(4) 



Q;(/i 



j(5) 



w 



w 



(30) 



where dij has been defined in 



and the LL running for a must be used. 
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Now we show the effect of the RG resummation on (^). Setting fiw = Mw, fJ'c = 
for the moment and expanding (|30D around to order a log (m^/M^) yields with (|29|): 



mm,)) 



167r^ 







1 + ^log 
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(0) / 



(31) 



Since J2 '''ijlij^ = 12(A^ — 1) / {2N), we have reproduced in (|3TD the large logarithm present 
in the coefficient Hll of (Bl). 



3.2 The Hamiltonian below the Scale fic 

At the scale /Xc we obtain C{nc) from the matching condition (^31) : 



-,(5) 



(32) 



The LL evolution of C(/i) below /ic can be obtained from the corresponding NLO expression 
in ( |2^ ) by setting Z^^^ to zero. 

Hence rji as defined in (|I]) or (^ reads in the LL approximation: 



(3) 



d 

(«(/^c)) - 

(a(/ic))^/^ 



/ a(/ife) 



^(4) 



3 / Q:(/ib) \ ^^^^^ ^ Q:(/ity 

2 1^^; 



\ 12/23 



_^1 ( aifib) 



-24/25 



2 V"(y"c 



a(/ib) j 

/ a(/xvy) 
V a{fib) 



'a(/i6 



-6/25 



,a(/ic 



\ ainb) 



-6/23 



-24/23' 



(33) 



In the preceding sketch of the LL calculation we have seen how the various scales /i^y, fJ^b 
and /ic enter the result of r]\^. H^, which is related to physical observables, must not depend 
on these scales, and we can vary them in (|33|) to judge the accuracy of the perturbative 
calculation. The residual scale dependence of ( p3| ) is sizeable and is substantially reduced 
in the NLO result. A comparative numerical analysis of the LL and NLO scale setting 
ambiguity will be postponed until section |, and at this place we just display the sources 
of the scale dependence analytically. 

The coupling constant grows with decreasing scale, is therefore most sensitive to 
the variation of the smallest scale /ic, and we exemplify the scale dependence by expanding 
the concerning parts of the hamiltonian around m* = mc{mc). They read 



7/}^«(/Xc)m^(/ic) 



{ml)' {a{m:))'- E 



aifib 



aim* 



a{jj, 



W) 



a{iJ,b) 
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a 
An 



log 



Iij 1 + 



/ n 



+ 0^2 log^ 



[ml 



[ml 



(34) 



Notice that the largest contribution to the coefficient of alog(/i^/m^) comes from the 

anomalous dimension of O The two-loop finite part of the NLO analysis will encounter 

a term which cancels the logarithm in (Rl) and will thereby diminish the scale ambiguity. 



4 The Next— to— Leading Order Analysis 

In this section we describe in detail our NLO calculation, in which the 0(a)-contribution 
to the transition amplitude (fig. |^) is improved by summing the terms log" (m^/M^), 
n = 0, 1, 2, . . . , to all orders in perturbation theory. 



4.1 The Hamiltonian between the Scales ji^ and jiw 

For the NLO factorization 

{H\,i)) = /xc</x</xiy, (35) 

we need the AS* = 1 Wilson coefficients C± to order a. At the matching point one has 



(36) 



with B± given in the NDR scheme 



B. = ±11^. (37) 



Hence in ( p5|) 



= l + %^(^log^ + 5. + 5,)+0(a^). (38) 



As outlined in section |2.1| the factorization ( pSf ) is so simple because of the absence of 
mixing between Oij and local operators. Let us briefly explain why they do not mix: The 
diagrams of fig. § and ^ contain two internal quarks with masses and mi, each of which 
is mu or mc- The elements of the anomalous dimension matrix responsible for the mixing 
can directly be read off from the divergent parts of these diagrams. Since the diagrams 
have dimension 2 and the anomalous dimension does not depend explicitly on in a mass 
independent renormalization scheme as MS, the singular part of the diagrams must be 
proportional to ml + mf by symmetry under ^ mi. Hence in the GIM-combination 
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{k, I) = (c, c) — 2(c, u) + {u, u) the divergent parts cancel and no mixing with local operators 
occurs. 

Since we had to calculate the two-loop contributions to {Oij) shown in fig. || for the 
matching at /ic described in section |4.2| , we could verify (|35|) and analyze the contributions 
of the individual diagrams to Cij: The diagrams D2, D4, Dg and D7 contain each one tree- 
level and one dressed AS = 1 four-quark operator as subdiagrams. All other diagrams 
can be divided by a three-particle cut through the gluon line and the u/c-quark lines into 
two subdiagrams, each of which contains a tree-level AS* = 1 operator with four (anti-) 
quarks and one gluon as external states. In ( ^5|) the latter match their corresponding 
counterparts of fig. |^ with the tree-level Wilson coefficient, thus these diagrams are also 
properly taken into account in (^). Hence there is no point in treating the diagrams Di, 
D3, D4 and D5 differently from the others as done by the authors of 0. In fact only D2, D4 

-part in (36). contains a dressed colour singlet AS" = 1 
-part of C± [O, and therefore matches 



and D7 contribute to the 0{a) 
operator, which does not contribute to the 0(a) 
the diagram Fg with the tree-level Wilson coefficient. Finally the term with log(/i^/M^) 
in ( |5B| ) originates solely from D7. We will see in a moment that this term will reduce the 
dependence of Cij on /iiy. 

The RG improved Wilson coefficient at /ic is immediately obtained from (^) and (38): 




(4) a{fiw 
^ij + 



a(/i6 



47r 



y(0 



+ B, + Bj 



(5) 



(39) 



We close this section by noting that in 
order a due to 



the LL scale ambiguity at fiw is removed to 



^(5) 

[a{iJ,w)] 



1 + 



(0) 



An 



log 



-,(5) 



4.2 The Hamiltonian below the Scale /ic 

To obtain the hamiltonian for < fic given in (^) we must first solve the matching 
condition between the three- and four-quark theory for C(/t): 

mf^c)) = ^xl E a,(/ic)(0.,(/ic)) = ^,xic{f^^){Q,M). (40) 

«j'=+,- 



As discussed in section this step requires the calculation of the finite parts of the 
two-loop diagrams shown in fig. ^. 
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4.2.1 Calculation of the Two— Loop Diagrams in Fig. |^ 

We now sketch some details of the two-loop calculation, which has been carried out inde- 
pendently by both of us. 

The diagrams of fig. ^ in arbitrary i?^-gauge for the gluon involve tensor integrals 
up to rank six. As explained in section |2.1| we could set the external momenta to zero, 
the masses rrik i of the internal quark lines connecting the four-quark operators were kept 
arbitrary and in the very end set equal to rric or to m„ = 0. Diagrams Di through D3 
are infrared singular and we have kept nonzero strange and down quark masses for their 
regularization. Alternatively one could use an off-shell external momentum or regulate 
both UV- and IR-singularities dimensionally. Needless to say that the same regularization 
must be used for all matrix elements in (^). It turned out to be convenient to calculate 
the integrals keeping the space-time dimension D arbitrary, which has resulted in quite 
compact expressions for the coefficients of the various tensors built out of combinations 
of metric tensors. These coefficients involve GauB' hypergeometric function 2F1, whose 
expansion in terms of e = (4 — L')/2 yields quite lengthy expressions containing the familiar 
dilogarithm function in the finite part. Di, D2 and D3 also contain a additive terms with 
the IR-regulators. The dilogarithms disappear for rrik = mi = rric or rrikj = mu — 0. Only 
a few number of the hypergeometric functions must be expanded explicitly, the others 
can be obtained via functional equations. As a byproduct we have obtained the results 
corresponding to a dimensional infrared regularization. The 1/e^-terms of the diagrams 
disappear in the GIM-combination. 

The D-dimensional Dirac algebra has been evaluated with the help of the package 
Tracer for the computer algebra system Mathematica. The results involve prod- 



ucts of up to five Dirac matrices on each fermion line. They have to be projected onto 
a set of Dirac structures which form a basis for D = 4. We have chosen the following 
convention: 

^^.jul&L^-f^'Yl^L = (16 - 4£) ® 7^L + ev. 
7p7W.7^^ ® 7^7 Vt = (256 - 224^) 7,L ® 7"L + ev. (41) 

where L = 1 — 75 and ev. stands for evanescent operators vanishing for D = 4. The infrared 
singular diagrams contain additional Dirac structures, which multiply only finite integrals. 

The same procedure has to be performed for the graphs with QCD counterterms for 
the inserted AS* = 1 operators. The diagrams of fig. ^ also contain divergent subloops 
involving both AS* = 1 operators, their divergent parts vanish in the GIM-combination 
like the divergence of fig. |], so that no counterterms for these sub divergences are needed. 

After summing the diagrams we indeed verify that the divergences disappear by GIM 



as stated in section |4.1| . We remark that the evanescent operators in ( ^iD can in general 



mix into physical operators WB. We have checked that this is not the case here. 



The results for the individual diagrams of fig. ^ and ^ can be found in the appendix. 
They represent 

(0,,(/i)) = (0,,(/i))(°) + ^(0.,(/i))«+0(a^) 
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with {Oij given in (H) and 



87r2 



(42) 



Let us analyze (P2| ) first by looking back at the factorization at = /iv^: The unphysical 
Dirac structures T and U come with the same coefficient functions hx and defined in 
(PUP as {H{nw))^^^ in (^. has been defined in (^), and with ^ij = 1 the factorization 
of the T- and [/-part at /i = /x^y becomes transparent in (pSf). We split a^i^'* in ( ^2]) into 
its physical part c^;^^ and those parts which depend on the infrared regulators or involve 
the gauge parameter: 



CLLif^) + nj ■ 



2 2 



2iV 



1 m:;m 
-Slog 



2N 
N-1 



log (mg//i^) -mglog(my/x^ 

9 2 

m„ — m j 



2A^ 



log 



(43) 



Clearly also in (^) the unphysical terms coincide with those of ([lOD , so that they correctly 
give no contribution in ( |55D to the Wilson coefficient. 
The desired physical parts in (|43| ) read: 



r,,3(l-iV)log^%) 



+ 



102 - 73A^ - 32iV2 + 3N^ 



-6 + 6N + N^ - 

12N ' 



r,_3(-l-iV)log^ 



+ 



12- 

34 - 39A^ + 8N^ - 3N^ 



8N 



-2 + AN- 3Ar2 ^ jys 
12iV ' 



r__ 3(-3 - A^) log 



-34 + 19iV + 12iV2 + 3N^ 
8N 



^2-6N + 5N^-N^ ,,,, 

+ vr^ . (44) 

12iV ^ ' 



With (^4|) one verifies ( pSj) as promised in section 



4.2.2 Determination of the Wilson Coefficient C{jj,) 

Next we want to use (|^) to solve (^0|) for C{^c)- To this end we also need {Qll), which 



is made up of the diagrams of fig. || and ^. One has flS 



(Qll(/x)) = (g^^(;x))W + ^[(g^^)(%(/i) + (f)(°)/.^ + (?7)(°)/.^ 



(45) 
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Here hx and hu have been defined in (0), and we extract the physical part c out of a(/x): 
'A^- 1 



a(/i) 



c + 



2N 



3 log- 



2 2 



Cf + 



N -1 
2N 



2-2 



log(m2//i2) _ ml log(my /i^ 



AT 



2iV 



1 m^m 
- log 



-3Cp-5 



2A^ 



(46) 



Now we are in a position to solve the matching condition (^) for C(/ic): 



C(/Xe) = m^(/ic) 



j(5) 



1 + 



a{Hb) - ajuc) ^(4) 



47r 



H : H : Oij{fiWj 



in 



An 



+ 4^ ^^Af^c 



(47) 



Here 



bijifJ'W) 



(0) 

log 



where i?± defined in 



originates from the matching at nw, while Tij and 



TijC 



(48) 



(49) 



have entered the result at the c-quark threshold /ic- All other quantities in (|47| ) are related 
to the RG evolution. 

Notice that again those parts of (Oij) and (Qll) which depend on the IR structure 
or the gauge parameter have not contributed to the Wilson coefficient. By working in an 
arbitrary i?^-gauge for the gluon we have explicitly proven that our final Wilson coefficient 
C{fi) is gauge independent. The use of small quark masses and as infrared regulators 
has enabled us to verify that C(yu) does not depend on the infrared structure. This would 
not have been so, if we had regulated both UV and IR singularities dimensionally, because 
the infrared structure would not have been explicit and the factorizations would have been 
trivial due to the absence of the operators T and If. 

C{fi) for fi < fic is readily obtained from (p^). The relation of C{fi) to rji is contained 
in (|1|), dH) and (||): 

C{fi) = r/i(/iH/,/ife,/ic)m^(/^c)&(/^), 
so that we can write down the final result now. 
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4.2.3 The Final Result for rji in Next— To— Leading Order 

The coefficient r/i in (|^) reads in next-to-leading order: 



(a(/ic)) 





a(/ic) (3) , aj^ib) - a{^c) ^(4) 
Att ^ Air ^ 



^.V H : bijinw) 



47r 



(50) 



The following table sums up in which equations the various quantities in (|50|) are 
defined: 













yif) 






m 


(20) 


(20) 


m 


m 


m 


m 


(|D,(|§,(|6|) 



Let us come back to the scale ambiguity at fic, which we had also discussed at the LL 
result by expanding ?7i (/ic)m^(yUc) around Hc = rn* = mdrnc) to order a log (/i^/(m*)^) 
in (0). In the NLO expression ([50|) we encounter via rij additional contributions to the 



coefficient of alog(/i^/(m*) ). They enter r^j in ( ^91) via c}^;^ given in (^41). Consider for 
example the coefficient of the logarithm in c^^. 1^ reads 



r__ 3 (-3 - A^) 



r 



7- 



(0) 



/ m 



and thereby cancels the corresponding coefficient in (P^. The same holds for the other 
terms in the sum so that in the NLO result (0) the scale ambiguity at /ic is removed to 
order a. 

We have expanded (0) further to estimate the remaining dependence on /ic. At order 
log (yU^/(m*)^) we find a large contribution of the ( — )-term, which originates in the 

fact that the anomalous dimension of O is negative and large in magnitude. In section 

^ we will discuss the scale ambiguities numerically. 

We close this section by mentioning another check of our result: If one switches off 

the running between and in (|50|) by expanding [a(/ivy)/a(/i;,)] etc. around a(/ivy), 

one must get the same result as by expanding 772 calculated in [0 in m^/M^ and setting 

rrit = rric. This is indeed the case. 



5 Numerical Results 

In this section we will point out the numerical impact of our NLO calculation on the value 
of Tji and the Kl - Ks mass difference. Special attention is paid to the dependence on the 
matching scales Hc and nw 
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The charm quark contribution to the AS = 2 Hamiltonian was given in (^. Replacing 
the S function therein by its lowest nonvanishing order in the expansion with respect to 
ml/My^ leads to 



r-<2 



167r2 



(51) 



with 



m 



*2 



Vlif'Wj /Xfe, /ic 



(52) 



rndfJ^c) needed in ( |5T| ) can be obtained from m* by using the mass evolution equation ([27|). 

H'^{fi) is directly related to physical observables and therefore must not depend on the 
scales fiw, fJ'b, fJ^c- For this reason we have introduced rj^ which multiplies only quantities 
which do not depend on these scales. t]1 is only relevant for the discussion of the dependence 
on /ic, since for fic = m* one has rji = r]*. 

Let us recall the expression for rji 



(3) 



[a 



j(5) 



w 



1 + 



a(/ic 



Air 



An 



I (3) a{nb) - ajuc) (4) 
-^^ + : ^ij 



Air 



+ 



An 



(53) 



and list its ingredients for the case = 3: 



(3) 



307 
"l62' 





++ 


+- = -+ 






12 

25 

12 

23 


6 

25 
6 

23 


24 
25 
24 
23 


7(4) 

7(5) 


6719 
1875 
5165 
1587 


419 
3750 
631 
3174 


7138 
1875 
4534 
1587 


Tij 


3 
2 


1 
2 


1 
2 




-91ogf -5-^ 


eiogf -f + 


_q incr !!i -u 35 _(_ Ze! 




41°g^ + f 




_8 loe - - 



20 



5.1 Numerical Results for 771 and ^uik 

This section lists the numerical values for rji and the - Ks mass difference Am^- The 
scale ambiguities will be viewed at in sect. |5l2 . 

For our numerical estimates we will use = A4 and the following list of input 



parameters [26 



Mw = 80 GeV, /i^ = m,, = 4.8 GeV, 

Bk = 0.7, rriK = 0.498 GeV, fx = 0.161 GeV, (Am)!^^ = 3.522 • 10"^^ GeV, 
ReAc = 0.215, Gp = 1.16639 ■ 10"^ GeV'l (54) 

In this section we will set the matching scales fic = f^t f^w = My/. Table || shows 
the LO and NLO result for rji for various values of rrtc and A^. 





0.150 


0.200 


0.250 


0.300 


0.350 




LO 


NLO 


LO 


NLO 


LO 


NLO 


LO 


NLO 


LO 


NLO 


1.25 


0.809 


0.885 


0.895 


1.007 


0.989 


1.154 


1.096 


1.334 


1.216 


1.562 


1.30 


0.797 


0.868 


0.877 


0.982 


0.965 


1.117 


1.064 


1.281 


1.175 


1.485 


1.35 


0.786 


0.854 


0.861 


0.960 


0.944 


1.085 


1.035 


1.235 


1.138 


1.419 


1.40 


0.775 


0.840 


0.847 


0.940 


0.924 


1.056 


1.010 


1.194 


1.105 


1.361 


1.45 


0.766 


0.828 


0.834 


0.922 


0.907 


1.030 


0.987 


1.157 


1.075 


1.310 


1.50 


0.757 


0.817 


0.822 


0.905 


0.890 


1.006 


0.966 


1.125 


1.048 


1.265 


1.55 


0.749 


0.806 


0.810 


0.890 


0.876 


0.985 


0.946 


1.095 


1.024 


1.225 



Table 1: r/i to LO and NLO for different values of A4 = Arg^ and = (both given in 
GeV). 

From this table we see, that the NLO correction is positive and sizable. For a typical 
value of A^ = 0.3 GeV it amounts to about 20 percent of the LO result. The dependence 
of rfi on Aj^ turns out to be stronger when NLO effects are included. This can be seen 
from fig. 0, where we have taken m* = 1.4 GeV. Larger values of Aj^ result in larger 
corrections and vice versa. The absolute value of rji decreases with increasing ttIc- 



Next, using (|5^), we calculate the contribution of to the Ki - Ks mass difference 



(Am)^ = ^niKflBKiReX^fnifi^l (55) 

with the result hsted in table 0. 

Since (Am)^ depends linearly on the Bk-, the result can easily be rescaled to other 
values of this parameter. 

Note that for m* = 1.4 GeV and A^jg = 0.3 GeV the contribution of to (Am)!^^ is 
about 64 percent. If we also take into account the top quark contribution in H^^^^ (|I|), 
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Figure 7: The dependence of rj^ on A4 = Aj^. 



which is about 6 percent, the complete short distance analysis reproduces about 70 percent 
of (Am)I^P. 

Using vacuum insertion = 1, the charm and top quark short distance contribution 
reproduces the full experimental mass difference, but [|TU], |TT], |T^ |T3[ favour a lower value 
for Bk- 



5.2 Scale Dependences of 771 

Let us now discuss the dependence of rj^ on the matching scale /ic- We will fix the values 
of mcirric) = m* = 1.4 GeV, /if, = 4.8 GeV and fiw = = 80 GeV. As can be seen from 
fig. H the scale dependence in the NLO calculation is reduced quite considerably compared 
to the LO result. 

If we take as a reference value = 0.3 GeV, then varying from 1.1 to 1.7 GeV 
amounts to a change of 63% in the LO case and 34% in the NLO calculation: 

r]l = 1.01 +Q2^ in LO , r/^ = 1.19 in NLO . 

For A^ = 0.2 GeV the corresponding variations are 45% in LO and 20% in NLO: 

rjl = 0.85 in LO , r/i = 0.94 ^Q-Jg in NLO . 
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Am 


Am 


Am 
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Am 


Am 


Am 


Am 


Am 


Am 


AmEXP 


AmEXP 


AmEXP 


AmEXP 


AmEXP 


1.25 


1.327 


0.377 


1.510 


0.429 


1.730 


0.491 


2.000 


0.568 


2.342 


0.665 


1.30 


1.408 


0.400 


1.593 


0.452 


1.812 


0.514 


2.078 


0.590 


2.409 


0.684 


1.35 


1.493 


0.424 


1.679 


0.477 


1.897 


0.539 


2.159 


0.613 


2.481 


0.705 


1.40 


1.580 


0.449 


1.768 


0.502 


1.986 


0.564 


2.245 


0.637 


2.560 


0.727 


1.45 


1.670 


0.474 


1.860 


0.528 


2.078 


0.590 


2.335 


0.663 


2.643 


0.751 


1.50 


1.763 


0.501 


1.955 


0.555 


2.173 


0.617 


2.428 


0.689 


2.731 


0.776 


1.55 


1.859 


0.528 


2.052 


0.583 


2.271 


0.645 


2.525 


0.717 


2.824 


0.802 



Table 2: (Am)^ to NLO in units of 10 GeV for different values of A4 = and 
/ic = rn* (both given in GeV). = 0.7. 



The variation of /ic between 1.1 and 1.7 GeV is quite conservative and may overestimate 
the 'theoretical error'. Note, that setting /ic ~ 1.25 GeV in the LO expression already 
reproduces the NLO result for fic = 1.4 GeV for all reasonable values of Ajj^. 

We will further analyse the remaining scale setting ambiguity by looking at the indi- 
vidual contributions coming from the different terms in the sum in ( |53D as we already did 
analytically in sec. |4.2.3| . 

In fig. ^one can easily see, that the (++) and (H — ) contributions behave essentially fiat 
as they should, but they contribute with opposite sign to rj^ and therefore almost cancel in 

the sum. On the other hand the ( ) part carries a large residual scale dependence and 

is the largest of the three contributions. 

Let us now continue with the discussion of residual scale dependences in rj* by varying 
the scale fiw We will fix fi^ = = 1.4 GeV, fib = 4.8 GeV, Mw = 80 GeV. Also here — 
as in the case of fic discussed above — the scale dependence is reduced considerably by the 
inclusion of NLO effects. Taking A^ = 0.3 GeV and varying from 60 to 100 GeV we 
find a change of 77^ by 7% in LO and 4% in NLO. This residual scale dependence therefore 
turns out to be small compared to the one found for fic- 

rj* turns out to be very stable when varying the fe-scale /if,. Even if we neglect completely 
the effects from an effective 5 quark theory by setting Hb = nw, therefore evolving from 
scale fiw down to fic in an effective 4 quark theory, the error turns out to be of the order 
of 1 percent. This provides us with a simple and useful approximation for the numerical 
evaluation of rji. 



i,j=+~ 




a(/ic) (3) ^ a{fiw) - ajf^c) y{4) 
Att Air ^ 



-Oij[fiw) 



+ ^rM] (56) 
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Figure 8: The dependence of rjl on the matching scale /ic, where the c quark is removed 
from the theory for = 0.1, 0.2, 0.3 GeV. m* — 1.4 GeV has been kept fixed. The other 
scales have been chosen to be — 4.8 GeV and nw — Mw — 80 GeV . 



6 Conclusions 

In the present work we have calculated the complete next-to-leading order short distance 
QCD corrections to the coefficient rji of the effective A^" = 2 hamiltonian describing the 
K^-K^ mixing. We have used dimensional regularization and the MS-scheme with anti- 
commuting 75. Large logarithms have been summed by factorizing the Feynman amplitudes 
and applying the rcnormalization group evolution to the Wilson coefficients. To this end 
the rcnormalization of bifocal operators had to be studied in detail. We have explicitly 
proven that all Wilson coefficients involved in the calculation depend neither on the gluon 
gauge nor on the infrared regulators. We have further checked that the dependence of the 
leading order result on the renormalization scales is removed at NLO to order a. 

The NLO short distance QCD corrections enhance r^i by 20 % with respect to the 
LO result. The short distance part of the K^-Ks mass difference is increased by the same 
amount, because it receives its dominant contribution from the first two quark generations. 
Finally the theoretical uncertainty due to the remaining dependence on the renormalization 
scale has been reduced by a factor of 2 compared to the LO result. 
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Figure 9: The dependence of the (++), (H — ) and ( ) contributions to t]^ on fi^ for 

^Ms = 0-3 GeV. The hue labeled by "all" is obtained by adding (++) + 2(+-) + ( ). 

Acknowledgements 

At first we would like to thank A. J. Buras for suggesting this project and for his per- 
manent encouragement during our work. We are grateful to him and G. Buchalla for 
carefully reading the manuscript and appreciate many helpful discussions with them, M. 
E. Lautenbacher and M. Misiak. 

A The Result for the Individual Diagrams 

In this appendix we collect the results for the individual diagrams depicted in fig. |^ and fig. 
^ including their QCD counterterms. The GIM mechanism has already been incorporated 
in the form (c, c) — {u,c) — (c, m) + {u,u), where (fc, /) denotes a diagram with internal k 
and / quarks. Then we get for the insertion of the bilocal operator O'^^ 

W = 0, (57) 
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where and denote colour factors listed in table |^. The gauge independent part 
and the one proportional to the gluon gauge parameter Df'^"^'^ read 



Dn 
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Df 



gauge 
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-3 Ri (m^,m^ 



+7^L (g) 7^L 
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Qgauge 
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+7^L (g) 7'^L 
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+7^L ® 7^L 
+7^L ® 7^^^ 
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-4 + ^-2 log(^) 



-2 - 3 log(^) 
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-2 + ^-log(4) 
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3 fj,^ 



(5J 



where the external spinors have been omitted. The mass singularities are contained in 

log (f)- log 



Ri ( 



2 ^ 



— m; 

o2 



Si (m2,m2 



^log(5j)-m^log 5 



(59) 



The expressions of diagrams generated from the generic ones in fig. ^ through a rotation by 
180 degrees or a left-right-refiection are obtained from eq. according to the following 
rules: 

Di^': 7^/2 ®7^L — > -f^L^^'^R, 



Dg : nis 



R^R 



To get the expressions for diagrams rotated by 90 degrees one has to apply a Fierz trans- 
formation. 
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1 


N^+2N^-2 


1 




1 


Af2-2 


1 


Ar3-2Af2+2 


2N 


2N 


2N 


2 


2N 


2 


2N 


2Af 


D5 


JV2-1 


N^+2N^-N-2 






JV2-1 


Ar2-i 


JV2-1 


N^-2N'^-N+2 


2N 


2N 


2N 


2 


2N 


2 


2N 


2N 


De 


N'^+N-l 
2N 


2N 


N^+N-l 
2N 


AT 
2 


N'^-N-1 
2N 


Af 
2 


2N 


Ar2-2 
2Af 


D7 


N-1 


N^+N-2 


N-l 


Af-1 


N+l 


Af+1 


N+l 


Ar2_7V-2 


2N 


2N 


2N 


2 


2N 


2 


2N 


2N 



Table 3: Color factors of the diagrams , ^ = 0, . . . , 7 
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